N92-I0760 


A 40™ DEGREE AND ORDER GRAVITATIONAL FIELD MODEL FOR MARS 

M.T. Zuber \ D.E. Smith\ F.J. Lerch \ R.S. Nerem \ G.fl. Patel 2 , and S.K. Frlcke 3 , 1 Laboratory for Terrestrial 
Dhyci^s, NASA/Qoddard Spscs Flight Csntsr, Qrssnbslt, MD 20771 ; ^ST Systsrns Corporation, Lanham, 
MD 20706; 3 RMS Technologies, Inc., Lanham, MD 20706. 


Understanding of the origin and evolution of major physiographic features on Mars, such as the 
hemispheric dichotomy and Tharsis rise, will require improved resolution of that planet’s gravitational and 
topographic fields. The highest resolution gravity model for Mars published to date [1 ] was derived from 
Doppler tracking data from the Mariner 9 and Viking 1 and 2 spacecraft, and is of 18™ degree and order. 
That field has a maximum spatial resolution of approximately 600 km, which is comparable to that of the 
best current topographic model [2]. The resolution of previous gravity models was limited not by data 
density, but rather by the computational resources available at the time. Because this restriction is no 
longer a limitation, we have re-anaiyzed the Viking and Mariner datasets and have derived a gravitational 
field complete to 40 th degree and order with a corresponding maximum spatial resolution of 300 km 
where the data permit. 

Derivation of the field was based on analysis of 235 orbital arcs consisting of 1 200 days of S-band 
Doppler tracking data from the Mariner 9 and Viking 1 and 2 spacecraft, collected by the Deep Space 
Network from 1971-1978. The data were processed using the GEODYN/SOLVE orbit determination 
programs. These programs, which have previously been used in the determination of a series of standard 
Earth gravitational models (the 'GEM* models [cf. 3]), have recently been adapted for analysis of planetary 
tracking data [4], 

The Martian gravitational potential at spacecraft altitude was represented in spherical harmonic 

form {5] 
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where r is the position vector of the spacecraft in areocentric coordinates, r is the radial distance from 

the center of mass of Mars to the spacecraft, <f> and X are the areocentric latitude and longitude of the 
spacecraft, r H is the mean radius of the reference ellipsoid of Mars, G is the gravitational constant, /Vf H is 
the mass of Mars, P lBj are the normalized associated Legendre functions of degree / and order m, C lm and 
S lm are the normalized spherical harmonic coefficients which were estimated from the tracking 
observations to define the gravitational model, and N is the maximum degree representing the size (or 
resolution) of the field. The gravitational force due to Mars which acts on the spacecraft corresponds to 
the gradient of the potential, V,. In our analysis, the origin of the field was taken to be the center of mass 
of Mars, which required that C 00 =1 and S 00 =C 10 =C 11 =S 11 =0. 

To determine the field, orbits were computed for each arc by estimating from the tracking data 
the initial position and velocity of the spacecraft, along with the atmospheric drag, solar radiation pressure, 
and Doppler tracking biases. After the arc solutions were iterated to convergence, information equations 
were created for each arc by evaluating the partial derivatives of the observations with respect to the arc 
parameters and gravity coefficients along each arc. The gravitational model was then found by adding 
together the information equations for each arc and solving the resulting linear system. The dominant error 
sources in the model are the uncertainties in the spacecraft orbits, which are affected by the tracking 
coverage as well as the assumed models of atmospheric drag and solar radiation pressure, and 
perturbations caused by spacecraft angular momentum desaturations. Details of the treatment of these 
parameters are discussed in [4]. As in previous analyses [e.g. 1], we imposed a priori constraints on the 
model based upon Kaula’s Rule [5] rescaled to Mars, which causes poorly observed (usually high degree 
and order) coefficients to tend toward zero, but has little effect on coefficients that are well sensed by the 
tracking data [3,6]. 

Error analyses based on the observation data demonstrate that the new field is characterized by 
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a significantly better rms orbital accuracy (=0.212 cm s' 1 ) than previous models ([1]=0.535 cm s' 1 ). In 
addition, the field shows a greater resolution of identifiable geological structures. Free air gravity 
anomalies calculated from the most recent iteration of the model, MGM-515, are plotted in Figure 1. As 
for previous models, the gravity anomalies correlate well with principal features of Martian topography, 
indicating either that topography is only partially compensated, or that deep-seated mantle processes are 
responsible for both the topography and gravity [7]. For the current model, gravity anomalies associated 
with essentially all major features, including volcanic shields, impact basins and the Valles Marineris, show 
considerably higher magnitudes than in previous models. A higher dynamic range of power (2100 m vs. 
1950 m for [1]) is also observed in the geoid representation of the field. Figure 1 also shows a possible 
trend along the hemispheric dichotomy of Mars in the longitude range 1 20° < A. <240°. Free air anomalies 
across the dichotomy boundary would be inconsistent with a simple model of isostatic compensation due 
to a change in crustal thickness across the boundary, such as suggested by [8], However, detailed 
modeling of individual arcs which cross the boundary in areas spatially removed from the influence of 
Tharsis and Elysium will be required to more accurately resolve the nature of the gravity signature of this 
feature. 
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Figure 1 . Free air gravity anomalies computed from Mars Gravity Model (MGM) 515 to 40x40. The contour 
interval is 50 mgals. 
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